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ABSTRACT: Fat quality is determined by the composition of fatty acids. Genetic relationships between this composition and single 
nucleotide polymorphisms (SNPs) in the stearoyl-CoA desaturase 1 (SCD1) gene were examined using 513 Korean native cattle. Single 
and epistatic effects of 7 SNP genetic variations were investigated, and the multifactor dimensionality reduction (MDR) method was 
used to investigate gene interactions in terms of oleic acid (C18:l), mono-unsaturated fatty acids (MUFAs) and marbling score (MS). 
The g. 6850+77 A>G and g. 14047 C>T SNP interactions were identified as the statistically optimal combination (C18:l, MUFAs and 
MS permutation p-values were 0.000, 0.000 and 0.001 respectively) of two-way gene interactions. The interaction effects of g. 6850+77 
A>G g.10213 T>C and g.14047 C>T reflected the highest training-balanced accuracy (63.76%, 64.70% and 61.85% respectively) and 
was better than the individual effects for C18:l, MUFAs and MS. In addition, the superior genotype groups were AATTCC, AGTTCC, 
GGTCCC, AGTCCT, GGCCCT and AGCCTT. These results suggest that the selected SNP combination of the SCD1 gene and superior 
genotype groups can provide useful inferences for the improvement of the fatty acid composition in Korean native cattle. (Key Words: 
Korean Native Cattle, Multifactor-Dimensionality Reduction (MDR), Single Nucleotide Polymorphism (SNP), Stearoyl-CoA 
Desaturase 1 (SCD1) Gene) 



INTRODUCTION 

Fat quality in beef market is determined by its fatty acid 
composition and fat tissue containing mono-unsaturated 
fatty acids (MUFAs), which have a lower melting point and 
thus contribute to the flavor and tenderness of beef 
(Miyazaki et al., 1979; Melton et al., 1982; Matsumoto et 
al., 2007). Oleic acid (CI 8:1) is one of the major MUFAs in 
beef fat and thought to be the source of the aroma in cooked 
beef (Melton et al., 1982; Mandell et al., 1998). For high- 
quality beef, it is important to improve its fatty acid 
composition, particularly its CI 8:1. In addition, the fatty 
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acid composition is a heritable trait with heritability ranging 
from 0.31 to 0.73 (Inoue et al., 2008). 

The stearoyl-CoA desaturase (SCD) gene, an enzyme in 
the chemical conversion of saturated fatty acids to 
unsaturated fatty acids, has been reported as a main factor 
for change of stearic acid (C18:0) and palmitic acid (C16:0) 
to C18:l on unsaturated fatty acid of w-9; and many studies 
have attempted to increase unsaturated fatty acids for the 
improved flavor of meat and health benefits (Kuchel et al., 
2004; Scollan et al., 2006; Alexander et al., 2007; Erkkila et 
al., 2008; Webb and O'Neill, 2008; Zhang et al., 2008). 
Recently, genetic relationships between the fatty acid 
composition of beef and multiple nucleotide sequence 
variants have been examined in encoding the SCD1 gene 
(Taniguchi et al., 2004). Three variants c.702 A<G, c.762 
T<C and c.878 C<T in the exon 5 region have been found 
to be associated with the composition of unsaturated fatty 
acids and MUFAs in the adipose tissue of Wagyu, Canadian 
Holstein, Jersey breeds, Fleckvieh bulls and Korean native 
cattle (p<0.05) (Kgwatalala et al., 2007; Milanesi et al., 
2008; Tsuji, 2008; Barton et al., 2010; Ohsaki et al., 2009; 
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Oh et al., 2011). SNPs in the UTR region of the SCD1 gene 
are significantly related to carcass and beef traits without an 
evaluation of fatty acids in Hanwoo (Shin et al., 2006). 
Based on the above discussion, the SCD gene, which is 
closely related to CI 8:1 as an important factor influencing 
the beef flavor, is considered to improve beef quality. This 
indicates that it is important not only to reveal the 
independent effect of each factor but also to analyze the 
combined effects. However, previous studies have paid little 
attention to single-factor effects and have not always 
considered independent and combined effects. 

The multifactor dimensionality reduction (MDR) 
method has been employed (Ritchie et al., 2001; Chung et 
al., 2006) as a technique for efficiently detecting multiple 
genes and their interaction effects. This method is designed 
to address high-dimensional data and uncover complex 
relationships without relying on specified models fitting 
interactions between multiple genes (Bastione et al., 2004). 
With the increased availability of SNP information from an 
increasing number of SNPs in livestock genome databases 
(www.animalgenome.org) as well as human and other 
mammalian genomes, functional SNPs represent a better 
option for the simultaneous identification of several SNPs 
based on high-throughput tools such as DNA chips 
(Barendse et al., 2007). 

In this study, the MDR method was used to detect 
superior SNP combinations influencing CI 8:1, MUFAs and 
marbling score (MS) in Korean native cattle. Given that the 
detection of genetic information on 15 SNPs in the intron, 
exon5 and 3 'UTR regions of the SCD1 gene are 
significantly associated with fatty acid composition and MS 
(Oh et al., 2011), the interaction effects of SNPs in the 
SCD1 gene on CI 8:1, MUFAs and MS were examined. 
Here multiple SNPs were used by employing the MDR 
method to test the main and interaction effects of multiple 
SNPs on the meat quality of Korean native cattle, and a 
number of SNP factors believed to affect CI 8:1, MUFAs 
and MS were considered. Finally, the superior genotype 
groups were explored based on interactions between 
candidate genes. 

MATERIALS AND METHODS 

Animals and the phenotype 

A total of 513 Korean native cattle (called Hanwoo in 
Korea) were bred in Gyeongsangbuk-Do, Korea. The cattle 
were fed according to the feeding program of each farmer 
(n = 14) and were the progeny of 18 sires. The pedigree 
records of these cattle were collected by the Charmpoom 
Hanwoo Research Institution. In general, they were weaned 
at 6 months of age, castrated at 6 months of age, fed 
growth-stage feed for 18 months, and given a high- 
concentration diet in the last 6 months. All were slaughtered 



at 941+72 d of age, and the marbling score was measured 
24 h after slaughter. First, the carcass was dissected at the 
last rib and the first lumber vertebra according to the 
Animal Product Grading System of Korea. The MS was 
measured in the left carcass cut across the vertebra between 
the last thoracic vertebra and the first lumbar vertebra. The 
marbling degree was scored from 1 to 9 (mean = 5.43) such 
that the higher the score, the more abundant the 
intramuscular fat. Genomic DNA was extracted from the 
longissimus muscle by using the LaboPass Tissue Mini kit 
(Cosmo Genetech, Seoul, Korea). 

Total lipids were extracted from approximately 500 mg 
of the longissimus muscle with chloroform/methanol (2/1, 
v/v) (Folch et al., 1957) and then methylated with sodium 
methylate (O'Keefe et al., 1968). They were filtered 
through a filter paper in a water bath (40°C). The filtrate 
was mixed with distilled water, from which a layer of 
methanol and water was removed. After the removal of 
chloroform and lipid layers by using nitrogen gas, the 
sample was treated with BF3-methanol (14%) and then 
subjected to transmethylation at 65°C. Fatty acid contents 
were analyzed using gas-chromatography (PerkinElmer, 
Inc., Waltham, MA, USA). 

SNP genotyping 

Fifteen SNPs of the SCD1 gene were selected based on 
the reference SNP of the bovine SCD1 gene (GenBank No. 
AY241932) in dbSNP. Genotypes at the 15 SNPs were 
subjected to a preliminary analysis. Primers for their 
amplification and extension were designed for single-base 
extension (SBE) (Vreeland et al., 2002) by using forward, 
reverse, and extension primer sequences (Supplementary 
Table 1). Primer extension reactions were performed using 
the SNaPshot ddNTP Primer Extension Kit (Applied 
Biosystems, Foster City, CA, USA). For the purification of 
the reaction of the primer extension, mixtures containing 
exonuclease 1 and shrimp alkaline phosphatase were added 
to the reaction mixtures. Samples were cultured at 37°C for 
1 h and then inactivated at 72°C for 15 min. The PCR 
product was mixed with the Genescan 120 LIZ standard and 
HiDi formamide (Applied Biosystems), followed by 
denaturation at 95 °C for 5 min. Electrophoresis was 
performed using the ABI PRISM 3130XL Genetic Analyzer, 
and then the electrophoresis product was assayed using 
GeneMapper v.4.0 (Applied Biosystems). 

Linear discriminant and statistical analyses 

For 45 extreme cattle, the significance of SNPs allele 
frequencies between the two groups (high and low) was 
analyzed using the chi-square test. Those SNP alleles for 
which there was a significant frequency difference were 
included in a linear discriminant analysis to develop a 
function classifying animals into the two groups (Tsuji et al., 
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Supplementary Table 1. The information of primer sequence about SNPs of SCD (Accession number : NM_AY241932) gene of BTA 
26 in Korean native cattle 



SNP Region 




Sequences 


Temperature 


Product 


g.6850+77 A>G Intron 


F 


ACTCTGGCCTGAGGAAAACA 


60 


348 




R 


CCGGAACTTAACCACAAGGA 








E 


GTGGCGGAATGGAGGCGATC 






g.8243-72 A>G 


F 


AGGCATCAGAATGGAACTGC 


60 


401 




R 


AGACACCCTCTCAGGGGAAT 








E 


AGCTGTGGGCTCTGGAAGAC 






g.8646+128 A>G 


F 


AGCTGGTGATGTTCCAGAGG 


60 


390 




R 


TTTGTGCCTCTCCTCGCTAT 








E 


AAATTACAATTTAAAACCCT 






g.l0153A>G Exon5 


F 


CAGAAAATTTCCTTGCCCATT 


60 


552 




R 


TGTTGCTTAACTTTCAAGGGTTT 








E 


ATCCTGCCCACACTCGTGCC 






g.10213 T>C 


F 


CAGAAAATTTCCTTGCCCATT 


60 


552 




R 


TGTTGCTTAACTTTCAAGGGTTT 








E 


TTTGCCACCTTATTCCGTTA 






g. 10329 C>T 


F 


CAGAAAATTTCCTTGCCCATT 


60 


552 




R 


TGTTGCTTAACTTTCAAGGGTTT 








E 


TCTGGTTTCCCTGGGAGCTG 






g. 13565 G>T 3TJTR 


F 


TTTTGTCCCTGTCACCTTCC 


60 


425 




R 


AGGGAGCATACTGGGTGATG 








E 


TCCAAAGCCTAGACAACCCT 






g. 13757 C>A 


F 


TTTTGTCCCTGTCACCTTCC 


60 


425 




R 


AGGGAGCATACTGGGTGATG 








E 


ACCAAACTGAAAAAACAGGT 






g. 14047 C>T 


F 


GTCTGCCCATCACCCAGTAT 


60 


397 




R 


TCACCACCATACCACCAATG 








E 


GAATTGTTTTCTTTCTTTCA 






g. 14578 A>G 


F 


AGAGGGGGCTTGAAGAAGAG 


60 


469 




R 


GCTTGCCTGTCCAGAAAAAG 








E 


GTTGTTGGGTGTCCTGATAT 






g. 15001 A>G 


F 


TGACAGAATGGTCAGGGTCA 


60 


328 




R 


TGACTCAGGGTAGGCAGCTT 








E 


ACTTGGCCTCCAGAGTCTCC 






g.15101 C>T 


F 


TGACAGAATGGTCAGGGTCA 


60 


328 




R 


TGACTCAGGGTAGGCAGCTT 








E 


GCATTCGCACTCATGAGGCA 






g.l5202G>A 


F 


TGACAGAATGGTCAGGGTCA 


60 


328 




R 


TGACTCAGGGTAGGCAGCTT 








E 


CAGAACCAGAGAGCATGCTG 






g. 15491 G>A 


F 


GGGAAGCACAAGAAGCAAAG 


60 


320 




R 


CAGTGGCAGAGCATTTGGTA 








E 


CAACACAGTCTAGGGCAGAC 






g. 16524 C>T 


F 


GGTAGTTTTAATTCCACAGGTCAT 


58 


332 




R 


TCCAAAATTATTCTTAGGTCACAA 








E 


ATGAAAAGTTCAAGAAAGCT 







1 Forward sequence. 2 Reverse sequence. 3 Extension sequence. 



2008). Only the influential fragments for discrimination (F 
>2.5) were included in the function. 

The Hardy-Weinberg equilibrium (HWE) was 
determined for each locus by comparing expected 
frequency with observed genotype frequencies by using the 



chi-squares test for the Korean cattle population. The 
relationships between phenotypes (C18:l, MUFAs and MS) 
and individual SNPs for the 513 samples were analyzed 
using the following mixed analysis of covariance 
(ANCOVA) model based on SPSS vl9.0 (SPSS Inc., 
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Chicago, IL, USA). 

yw = M+Si+Pj+Gk+fiA t e+e<ik 

where y ijU is the phenotype, jU is the overall mean, 5, 

is the random effect of sire i with an assumption of 
independent and identical normal distribution, Pj is the 
fixed effect of calving place j (14 classes), Gu is the fixed 
effect of genotype k, ft is the regression coefficient, age is a 
covariate for age in days at slaughter, egu is the random 
residual assumed to have independent and identical normal 
distribution. 

MDR analysis 

The MDR method, a non-parametric and model-free 
technique, was initially implemented in case-control studies 
(Hahn et al., 2003). For its application to continuous data, 
the classification and regression tree (CART) algorithm has 
been developed and combined with the MDR method 
(Paolo, 2003). The expanded MDR method involves 
classification into two groups using a regression tree, i.e., 
high and low groups for phenotypes, and the group's 
impurity can be defined as: 



i(g) = 



Where y g is the mean for the node (group) and y ffl is 

the observation of the j th individual in a total of N g 
individuals in the group (g = high or low). Each individual 
is assigned to a cell (e.g., if there are two SNPs with three 
genotypes per SNP, then there are 9 possible cells), and then 
each cell can be defined as a high or low group. The 
expected numbers of high and low groups are: 



N hish = £ I hish ii) ■ N t • and N low = X I„Ji) ■ N, 



Respectively, where n is the number of total cells and 
I g (i) is an indicator function in which 



1 i(cell) s high group f and 

0 o.w 

1 i(cell) e low group 
0 o.w 



Then the average squared error (ASE) is: 

ASE = Shigh/Nhigh+SioJNioy, 



where e 



high ' 



5>, 



is the sum 



of squared errors for 

N, 



the 



high group, 



is that for low group, 



N, 



and yij is the phenotype of the j' h individual in the i' h cell. 

In this study, the procedures for the expanded MDR 
method were as follows: 

Step 1: Data were randomly divided into 10 equal parts 
(one testing set and nine training sets for cross-validation 
purpose). 

Step 2: A set of n SNPs was selected from the pool of all 
SNPs. 

Step 3: Based on the observed level of each n, steers 
were partitioned into classes, referred to as cells. If n = 2, 
then a set of two SNPs was selected, and, because a SNP 
had three genotypes, there were 3 2 = 9 possible cells. 
Phenotypic means were calculated within each cell. 

Step 4: The impurity function in the CART algorithm 
used the variance impurity such that the group with the 
higher average value was labeled as high and the rest was 
labeled as low. 

Step 5: The expanded MDR model with the smallest 
ASE was chosen among all two-factor combinations. 

Step 6: For the evaluation of the predictive ability of the 
model, the predicted ASE (P_ASE) was estimated using a 
10-fold cross-validation method. These six steps were 
repeated for each possible combination of given n. The 
model with the minimum predicted ASE was selected as the 
best model, but for the selected model, statistical 
significance was not determined by the predicted ASE. 
Therefore, permutation tests were performed to determine 
the empirical significance thresholds by applying the same 
MDR method (Good, 1994). Before the MDR 
implementation, phenotypes were adjusted for 
contemporary and, sire/steer's age effects by using residuals 
obtained after fitting the general linear model without SNP 
effects. 

RESULTS AND DISCUSSION 

The SCD is an endoplasmic reticulum (ER) enzyme that 
catalyzes the biosynthesis of MUFAs from saturated fatty 
acids that are either synthesized de novo or derived from the 
diet (Smith et al., 2006). Variations in SCD enzyme activity 
in mammals are likely to affect a variety of key 
physiological variables, including cellular differentiation, 
insulin sensitivity, metabolic rates, adiposity, atherosclerosis, 
cancer, and obesity (Paton and Ntambi, 2009). Two types of 
SCD gene isoforms have been characterized in cattle, 
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namely SCD1, expressed mainly in adipose tissue, and 
SCD5, expressed mainly in the brain (Lengi and Corl, 
2007). The expression of SCD1 in bovine adipose tissue is 
regulated by numerous factors. The SCD1 mRNA level in 
muscle tissue has been reported to increase after weaning 
until 12 mon of age (Lee et al., 2005) but also during the 
late fattening stages (Kwon et al., 2009), and its 
development was dependent on the breed and diet (Chung 
et al., 2007). Breed differences in the SCD1 mRNA level 
have been found Japanese Black and Holstein steers 
(Taniguchi et al., 2004). Based on these reports, we 
collected 15 SNPs suggested by dbSNP of NCBI to identify 
those SNPs related to fatty acids in the SCD1 gene (Oh et 
al., 2011). These SNPs were located at 3 SNPs in the intron 
region, 3 SNPs in the exon 5 region and 9 SNPs in the 
3'UTR region. By SBE method for verification of 15 SNPs 
in 45 extreme steers for the high and low MUFAs group of 
Korean native cattle, 15 polymorphic SNPs were identified 
(Table 1). And they were excluded in the following genetic 
association analysis. The other loci had their minor allele 



frequency ranging from 0.233 to 0.500, and their genotypes 
did not deviate from the Hardy-Weinberg equilibrium 
(p>0.05, Table 1). 

Table 2 shows the efficiency of discrimination for the 15 
polymorphic SNPs in the high and low MUFAs groups (n = 
45). For the classification of 45 extreme animals into these 
two groups, the discriminant coefficients were calculated 
for the 15 polymorphic SNPs by using information on SNP 
genotypes the high- and low-MUFAs groups. The 
discriminant coefficients for 7 of these 15 SNPs (g. 6850+77 
A>G, g.8646+128 A>G, g.10153 A>G, g.10213 T>C, 
g. 10329 C>T, g. 14047 C>T and g. 14578 A>G) were -1.786, 
-0.500, -1.267, -0.424, -0.954, -0.966, and -0.944, 
respectively. These 7 SNPs showed significant relationships 
with MUFAs (p<0.05, F>2.5). This is equivalent to a 
discriminative error of 2.1% and a high correlation ratio of 
0.91. These results indicate that the discriminant function 
successfully separated extreme cattle (Tsuji et al., 2008). 
Therefore, these 7 SNPs were selected for large scale 
genotyping in the Korean native cattle (n = 513). 



Table 1. Genotype and frequency of 15 polymorphic SNPs of the SCD gene in Korean native cattle 



SNP 



Region 



Genotype (No. of animals) 



Frequency 



H 



MAF 



HWE 



g.6850+77 A>G Intron 


AA(2) 


GA(17) 


GG (26) 


N 4 (45) 


0.358 


0.244 


0.708 




0.044 


0.378 


0.578 


1 








g.8243-72 A>G 


AA(3) 


GA (22) 


GG (20) 


N(45) 


0.429 


0.233 


0.345 




0.067 


0.489 


0.444 


1 








g.8646+128 A>G 


AA(9) 


GA(21) 


GG(15) 


N(45) 


0.491 


0.311 


0.738 




0.200 


0.467 


0.333 


1 








g.l0153A>G Exon 5 


AA(10) 


GA(19) 


GG(16) 


N(45) 


0.491 


0.433 


0.345 




0.222 


0.422 


0.356 


1 








g. 10213 T>C 


CC(13) 


CT (22) 


TT(10) 


N(45) 


0.498 


0.448 


0.905 




0.289 


0.489 


0.222 


1 








g. 10329 C>T 


CC (10) 


CT (20) 


TT(15) 


N(45) 


0.494 


0.467 


0.502 




0.222 


0.444 


0.333 


1 








g.l3565 G>T 3'UTR 


GG (8) 


GT(21) 


TT(16) 


N(45) 


0.484 


0.444 


0.808 




0.178 


0.467 


0.356 


1 








g. 13757 C>A 


AA(12) 


GA(22) 


GG(ll) 


N(45) 


0.499 


0.489 


0.884 




0.267 


0.489 


0.244 


1 








g. 14047 C>T 


CC (10) 


CT (23) 


TT(12) 


N(45) 


0.499 


0.478 


0.871 




0.222 


0.511 


0.267 


1 








g. 14578 A>G 


AA(ll) 


GA(20) 


GG (14) 


N(45) 


0.498 


0.467 


0.472 




0.244 


0.444 


0.311 


1 








g. 15001 A>G 


AA(10) 


GA(20) 


GG(15) 


N(45) 


0.494 


0.444 


0.502 




0.222 


0.444 


0.333 


1 








g.15101 C>T 


CC (10) 


CT (23) 


TT(12) 


N(45) 


0.499 


0.478 


0.871 




0.222 


0.511 


0.267 


1 








g.l5202G>A 


AA(ll) 


GA(23) 


GG(ll) 


N(45) 


0.500 


0.500 


0.881 




0.244 


0.511 


0.244 


1 








g. 15491 G>A 


AA(13) 


GA(23) 


GG (9) 


N(45) 


0.496 


0.456 


0.839 




0.289 


0.511 


0.200 


1 








g. 16524 C>T 


CC (8) 


CT (26) 


TT(ll) 


N(45) 


0.498 


0.467 


0.280 




0.178 


0.578 


0.244 


1 









Heterozygosity. 2 Minor allele frequency. 3 Hardy-Weinberg equilibrium. 4 Total number. 
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Table 2. Efficiency of discrimination for 15 ploymorphic SNPs in high and low MUFAs groups of Korean native cattle (n = 45) 



SNP 



High MUFAs group 



Low MUFAs group 



observed genotype 
No. of animals (Frequency) 



observed genotype 
No. of animals (Frequency) 



Discriminant 
coefficient 



F-value p-value 



g.6850+77 A>G 


AA 


GA 


GG 


AA 


GA 


GG 


-1.786 


5.041 


0.011 




2(1.00) 


12(0.71) 


8(0.31) 


0(0.00) 


5(0.29) 


18(0.69) 








g.8243-72 A>G 


AA 


GA 


GG 


AA 


GA 


GG 


-0.305 


0.241 


0.787 




2(0.67) 


11(0.5) 


9(0.45) 


1(0.33) 


11(0.5) 


11(0.55) 








g.8646+128 A>G 


AA 


GA 


GG 


AA 


GA 


GG 


-0.500 


3.041 


0.058 




4(0.44) 


14(0.67) 


4(0.27) 


5(0.56) 


7(0.33) 


11(0.73) 








g.10153 A>G 


AA 


GA 


GG 


AA 


GA 


GG 


-1.267 


5.358 


0.008 




7(0.70) 


12(0.63) 


3(0.19) 


3(0.30) 


7(0.37) 


13(0.81) 








g.l0213T>C 


TT 


CT 


CC 


TT 


CT 


TT 


-0.424 


3.907 


0.028 




4(0.40) 


15(0.68) 


3(0.23) 


6(0.60) 


7(0.32) 


10(0.77) 








g. 10329 C>T 


cc 


CT 


TT 


CC 


CT 


TT 


-0.954 


4.255 


0.021 




6(0.60) 


13(0.65) 


3(0.20) 


4(0.40) 


7(0.35) 


12(0.80) 








g. 13565 G>T 


GG 


GT 


TT 


GG 


GT 


TT 


0.185 


0.248 


0.781 




3(0.38) 


11(0.52) 


8(0.50) 


5(0.62) 


10(0.48) 


8(0.50) 








g. 13757 C>A 


CC 


CA 


AA 


CC 


CA 


AA 


-0.811 


1.919 


0.159 




8(0.73) 


10(0.46) 


4(0.33) 


3(0.27) 


12(0.54) 


8(0.67) 








g. 14047 C>T 


CC 


CT 


TT 


CC 


CT 


TT 


-0.966 


3.739 


0.032 




6(0.60) 


14(0.61) 


2(0.17) 


4(0.40) 


9(0.39) 


10(0.83) 








g. 14578 A>G 


AA 


GA 


GG 


AA 


GA 


GG 


-0.944 


3.337 


0.045 




7(0.64) 


12(0.60) 


3(0.21) 


4(0.36) 


8(0.40) 


11(0.79) 








g. 15001 A>G 


AA 


GA 


GG 


AA 


GA 


GG 


-0.069 


0.021 


0.979 




5(0.50) 


10(0.50) 


7(0.47) 


5(0.50) 


10(0.50) 


8(0.53) 








g.15101 C>T 


CC 


CT 


TT 


CC 


CT 


TT 


0.179 


0.199 


0.821 




4(0.40) 


12(0.52) 


6(0.50) 


6(0.60) 


11(0.48) 


6(0.50) 








g. 15202 G>A 


GG 


GA 


AA 


GG 


GA 


AA 


0.000 


0.095 


0.909 




5(0.46) 


12(0.52) 


5(0.46) 


6(0.54) 


11(0.48) 


6(0.54) 








g. 15491 G>A 


GG 


GA 


AA 


GG 


GA 


AA 


0.185 


0.098 


0.907 




4(0.44) 


11(0.48) 


7(0.54) 


5(0.56) 


12(0.52) 


6(0.46) 








g. 16524 C>T 


CC 


CT 


TT 


CC 


CT 


TT 


-0.302 


0.343 


0.712 




5(0.63) 


12(0.46) 


5(0.46) 


3(0.37) 


14(0.54) 


6(0.54) 









The MDR method, a novel technique, was applied to 
examine the interactions between the 7 polymorphic SNP 
genotypes in complex traits (C18:l, MUFAs and MS). 
These were applied to the MDR method, and superior two- 
and three-way SNP combinations were used for each trait. 
Here p<0.05 was used to determine statistical significance 
for the critical value in the t-test and permutation test. For 
the investigation of two- and three-way SNP gene 
interactions, the top 5 combinations were listed for each 
trait (Tables 3 and 6). The ranking was determined based on 
the training-balanced accuracy of results for each trait. 

Table 3 shows that the g. 10213 T>C and g. 14047 C>T 
combination had the highest score for both training- 
balanced accuracy (60.16%) and testing -balanced accuracy 
(57.42%) in the two-way interaction model of CI 8:1. The 
g. 10329 C>T and g. 14047 C>T combination represented 
the best two-way interaction model of MUFAs. In case of 
MS, the g.10153 A>G and g.10213 T>C combination 



represented the highest training-balanced accuracy 
(56.78%). 

The best interaction models of CI 8:1, MUFAs and MS 
were not the same, and therefore 4 superior SNP 
combinations (g.10213 T>C and g.14047 C>T; g.14047 
C>T and g. 14578 A>G; g. 10329 C>T and g.14047 C>T; 
and g.6850+77 A>G and g.14047 C>T) were calculated and 
compared (Table 4). Then, the g.6850+77 A>G and g.14047 
C>T SNP interaction was selected as the statistically best 
combination (C18:l, MUFAs and MS permutation t-values 
were 7.091, 7.017 and 3.293 respectively) of two-way 
interactions. 

The single and combination effects of g.6850+77 A>G 
and g.14047 C>T were examined (Table 5). As shown in the 
table, the combination effects of g.6850+77 A>G and 
g.14047 C>T on C18:l, MUFAs and MS traits were better 
than their individual effects. In particular, g.6850+77 A>G 
and g.14047 C>T SNPs were not significant (with 
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Table 3. Significant SNP combinations of top 5 two-way gene interactions from 3 traits based on MDR method in Korean native cattle 







C18:! 1 




MUFAs 2 




MS 3 






Training- 


Testing- 


olNr 


Training- 


Testing- 


oINr 


Training- 


Testing- 


balanced 


balanced 


balanced 


balanced 


balanced 


balanced 






accuracy 


accuracy 




accuracy 


accuracy 




accuracy 


accuracy 


l 


g .10213T>C, 


0.6016 


0.5742 


g. 10329 OT, 


0.6140 


0.5585 


g.l0153A>G 


0.5678 


0.4839 




g. 14047 OT 






g. 14047 OT 






g. 10213 T>C 






2 


g. 14047 OT, 


0.5997 


0.5625 


g. 14047 OT, 


0.6131 


0.5644 


g.6850+77 A>G 


0.5642 


0.4664 




g. 14578 A>G 






g. 14578 A>G 






g. 10329 OT 






3 


g. 10329 C>T, 


0.5386 


0.5640 


g.6850+77 A>G, 


0.6109 


0.5387 


g.6850+77 A>G 


0.5634 


0.4675 




g. 14047 OT 






g. 14047 OT 






g,14578A>G 






4 


g.l0153A>Q 


0.5968 


0.5683 


g. 10213 T>C, 


0.6056 


0.5821 


g.6850+77 A>G 


0.5608 


0.4827 




g. 14047 OT 






g. 14047 OT 






g. 14047 OT 






5 


g.6850+77A>G 


0.5961 


0.5464 


g. 8646+1 28 A>G, 


0.6013 


0.5762 


g.8646+128A>G, 


0.5575 


0.4927 




g. 14047 OT 






g. 14047 OT 






g,10213T>C 







1 Oleic acid. 2 Monounsaturated fatty acids. 3 Marbling score. 4 SNP combination. 



permutation p values of 0.060 and 0.192 respectively) in 
MS. However, the g.6850+77 A>G and g. 14047 C>T 
combination was highly significant with p-value 0.001. In 
terms of superior genotypes CCAA, CCAG CCGG and 
CTAG were obtained from C18:l (45.24+0.18); CCAA, 
CCAG CTAA and CTAG from MUFAs (54.70±0.23) and 



CCAA, CCAG CCGG and CTAG from MS (5.76±0.13). 
Therefore, the CCAA, CCAG and CTAG genotype of the 
g.6850+77 A>G and g. 14047 C>T combination had the 
greatest interaction effects on all 3 traits. Because all these 
combinations were significant, the best combination was 
determined based on the t-value. 



Table 4. The selected of 4 combinations based on two-way gene interactions for 3 traits in Korean native cattle 



SNP 1 


Traits 2 


Genotype 


N 3 


LSmean±SE 


t-value 


Permutation 

P 4 


g.l0213T>C, g.14047 OT 


C18:l 


CCTT, CCTC, CTTC, CTCC 


315 


44.92±0.14 


6.921*** 


0.000 






Others 


198 


43.32±0.18 








MUFA 


CCTT, CCTC, CTTC, CTCC 


315 


54.17±0.16 


6.691*** 


0.000 






Others 


198 


52.4410.20 








MS 


CCTT, CCTC, CTTT, CTTC, CTCC, 


397 


5.5610.10 


2.879** 


0.005 






TTTT, TTTC 














Others 


116 


4.9710.17 






g. 14047 OT,g. 14578 A>G 


C18:l 


AACC, AGCC, AGCT 


294 


44.9710.15 


6.941*** 


0.000 






Others 


219 


43.4010.17 








MUFA 


AACC, AACT, AGCC, AGCT 


300 


54.2410.16 


6.952*** 


0.000 






Others 


213 


52.4710.19 








MS 


AACC, AGCC, AGCT, AGTT, GGCC 


326 


5.5910.11 


2.560* 


0.014 






Others 


187 


5.1410.14 






g. 10329 OT, g. 14047 C>T 


C18:l 


CCCC, CCCT, CTCT, CTTT 


318 


44.9110.14 


6.888*** 


0.000 






Others 


195 


43.3110.18 








MUFA 


CCCC, CCCT, CTCT, CTTT 


318 


54.1810.16 


6.828*** 


0.000 






Others 


195 


52.4110.20 








MS 


CCCC, CCCT, CTCT, CTTT, TTCC, 


373 


5.6010.10 


3.290** 


0.002 






TTCT 














Others 


140 


4.9710.16 






g.6850+77 A>G g. 14047 OT 


C18:l 


CCAA, CCAG CCGG CTAG 


214 


45.2410.18 


7 091*** 


0.000 






Others 


299 


43.6310.14 








MUFA 


CCAA, CCAG CTAA, CTAG 


181 


54.7010.23 


7 017*** 


0.000 






Others 


332 


52.8610.15 








MS 


CCAA, CCAG CCGG CTAG, TTAA, 


214 


5.7610.13 


3.293** 


0.001 






TTAG 














Others 


299 


5.1910.11 







SNP combination. C18:l = Oleic acid, MUFA= Monounsaturated fatty acids, MS = Marbling s 
3 Number of animals. 4 Permutation p value. * p<0.05. ** p<0.01. *** p<0.001. 
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Table 5. Sin; 


jle and combinations effect of 2 SNPs for 3 traits in Korean native cattle 








Traits 1 


SNP 


frPn ntvnp 


N 2 


LSmean±SE 


t-value 


p 3 


C18:l 


p 68^0+77 A>G 


AA, AG 


220 


45.03±0.20 


5.380*** 


0.000 






GG 


293 


43.75±0.14 








p 1 4047 C>T 


CC, CT 


332 


44.86±0.14 


6.688*** 


0.000 






TT 


181 


43.28±0.19 








kJll 1 Ulll tX 11 \J 11 


ccaa cc ag rrnn TTAG 


214 


45.24±0.18 


7.091*** 


0.000 






w Liici a 


299 


43.63±0.14 






MUFAs 


p 68^0+77 A>G 


AA, AG 


220 


54.28+0.22 


5.136*** 


0.000 






GG 


293 


52.92±0.15 








p 1 4047 G^.T 

g.l^tU^T/ V J 1 




%V) 


54.09±0.16 


6.268*** 


0.000 






TT 
1 1 


1 S 1 
lot 


52.43±0.21 








Combination 


CCAA, CCAG CTAA, CTAG 


181 


54.70±0.23 


7 017*** 


0.000 






WLI1CIS 




52.86±0.15 






MS 


g.6850+77A>G 


AA, AG 


220 


5.29±0.13 


1.802 


0.060 






GG 


293 


5.60±0.11 








g. 14047 C>T 


CC, CT 


332 


5.24±0.11 


1.348 


0.192 






TT 


181 


5.51±0.15 








Combination 


CCAA, CCAG CCGG CTAG, 


214 


5.76±0.13 


3.293** 


0.001 






TTAA, TTAG 














Others 


299 


5.19+0.11 







1 C18:l = Oleic acid, MUFA = Monounsaturated fatty acids, MS = Marbling score. 

2 Number of animals. 3 Permutation p value. * p<0.05. ** p<0.01. *** p<0.001. 



For the three-way interaction model, the g. 6850+77 
A>G, g. 10213 T>C and g.14047 C>T combination showed 
the highest training-balanced accuracy (63.76%, 64.70% 
and 61.85% respectively) (Table 6). Noteworthy is that 
g.6850+77 A>G and g. 14047 OT were also observed in 
the two-way interaction model. As shown in Table 7, the 
individual and interaction effects of g.6850+77 A>G, 
g.10213 T>C and g. 14047 OT were examined and 



compared for CI 8:1, MUFAs and MS. The superior 
genotype groups were 45.45+0.13 for CI 8:1, 54.73+0.22 
for MUFAs and 5.95+0.12 for MS and significant (all 
permutation p-values were 0.000) (Table 7). By contrast, 
the average for the htl*ht2 haplotype group (g. 10153 A>G, 
g.10213 T>C and g.10329 C>T) was 44.92 for C18:l, as 
described in Oh et al. (2011); However, in this study, the 
average for genotype of g.6850+77 A>G (AA, AG), 



Table 6. Significant SNP combinations of top 5 three-way gene interactions from 3 traits based on the MDR method in Korean native 
cattle 







C18:l' 




MUFAs 2 




MS 3 


Rank 


SNP 4 


Training- 


Testing- 


SNP 


Training- 


Testing- 


SNP 


Training- 


Testing- 


balanced 


balanced 


balanced 


balanced 


balanced 


balanced 






accuracy 


accuracy 




accuracy 


accuracy 




accuracy 


accuracy 


1 


g.6850+77 A>G 
g.10213 T>C, 
g. 14047 C>T 


0.6376 


0.5851 


g.6850+77 A>G 
g.10213 T>C, 
g. 14047 OT 


0.6470 


0.5789 


g.6850+77 A>G 
g.10213 T>C, 
g. 14047 OT 


0.6185 


0.5900 


2 


g.8646+128A>G 
g.l0153A>G 
g. 14047 C>T 


0.6351 


0.5860 


g.8646+128 A>G 
g. 10153 A>G 
g. 14047 OT 


0.6432 


0.5758 


g.6850+77 A>G 
g.8646+128 A>G 
g.10213 T>C 


0.6092 


0.5722 


3 


g.6850+77 A>G 
g.l0153A>G 
g. 14047 OT 


0.6314 


0.5927 


g.6850+77 A>G 
g. 10329 OT, 
g. 14047 OT 


0.6401 


0.5619 


g.6850+77 A>G 
g.10213 T>C, 
g. 13757 OA 


0.6072 


0.5608 


4 


g.6850+77 A>G 
g. 10329 C>T, 
g. 14047 OT 


0.6287 


0.5923 


g.8646+128 A>G 
g. 10329 OT, 
g. 14047 OT 


0.6339 


0.6117 


g.6850+77 A>G 
g.10213 T>C, 
g. 14578 A>G 


0.6057 


0.5737 


5 


g.6850+77 A>G 
g. 14047 OT, 
g. 14578 A>G 


0.6233 


0.5608 


g. 13757 OA, 
g. 14047 OT, 
g. 14578 A>G 


0.6310 


0.6263 


g.6850+77 A>G 
g. 10153 A>G 
g. 10329 OT 


0.5979 


0.5586 



1 Oleic acid. 2 Monounsaturated fatty acids. 3 Marbling score. 4 SNP combination. 
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Table 7. Single and combinations effect of 3 SNPs for 3 traits in Korean native cattle 



Traits 1 SNP Genotype N 2 LSmean±SE t-value p 3 



C18:l g.6850+77A>G 


AA, AG 


220 


45 03+0 20 


5.380*** 


0.000 




GG 


293 


43.75+0.14 






g.10213 T>C 


TT 


112 


44.91 ±0.27 


2.743** 


0.013 




TC, CC 


401 


44 1 "3+0 1 3 






g. 14047 C>T 


CC, CT 


332 


44 86+0 1 4 


6.688*** 


0.000 




TT 


181 


4^ 78+0 19 






Combination 


AATTCC, AGTTCC, AATCCC, 


189 


4^ 4S+0 1 3 


7.484*** 


0.000 




AGTCCC, GGTCCC, AGTCCT, 












GGCCCT, AGCCTT 












Others 


324 


43.63±0.20 






MUFAs g.6850+77 A>G 


AA, AG 


220 


54 78+0 77 


5.136*** 


0.000 




GG 


293 


S? Q9+0 1 S 






p 1 0? 13 T>r 


TT TC 


393 


« 78+0 1 s 

JJ. / OlLU. 1 J 


3.796*** 


0.003 




CC 


120 


JZ.OZltU.Z / 






a 1 4047 C>T 


CC, CT 


332 


r»Q+n i ^ 

JH- .VyAX). ID 


6.268*** 


0.000 




TT 


181 








n m r> i n a 1 1 n n 

V— ^ k_Jll 1 L/lll Cl LI W 11 


AATTCC AGTTCC AATCCC 


204 


7^+0 11 


7.768*** 


0.000 




AGTCCC, GGTCCC, GGTTCT, 












AATCCT, AGTCCT, GGCCCT, 












AGCCTT 












Others 


309 


52.70±0.15 






MS g.6850+77 A>G 


AA, AG 


220 


s fso+o 1 ^ 


1.802 


0.060 




GG 


293 


S 7Q+H 1 1 
J.Z71U. 1 1 






a i o? 1 ^ tsp 


TT TC 


^Q^ 

jyj 


j.jo±U.lU 


7 878** 


0 006 




CC 


120 


4.99+0.17 






g. 14047 C>T 


CC, CT 


332 


5.51+0.11 


1.348 


0.192 




TT 


181 


5.27±0.15 






Combination 


AATTCC, AGTTCC, GGTCCC, 


237 


5.95±0.12 


5.884*** 


0.000 




AATTCT, GGTTCT, AGTCCT, 












AACCCT, AGCCCT, GGCCCT, 












AATTTT, GGTTTT, AATCTT, 












GGTCTT, AACCTT, AGCCTT 












Others 


276 


4.97+0.11 







1 C18:l = Oleic acid, MUFA= Monounsaturated fatty acids, MS = Marbling score. 

2 Number of animals. 3 Permutation p value. * p<0.05. ** p<0.01. *** p<0.001. 



g.10213 T>C (TT, TC) and g.14047 C>T (CC, CT) SNP 
combinations was 45.45, better than 0.53 for the htl*ht2 
haplotype group. The combinations effect of SNPs was 
greater than the single effect for CI 8:1, MUFAs and MS 
(Table 7). As in the case of the two-way interaction effect, 
the individual effect of g.6850+77 A>G and g. 14047 C>T 
on the MS were not significant. In addition, g.6850+77 A>G, 
g.10213 T>C and g. 14047 C>T had a significant three-way 
interaction effect, and AATTCC, AGTTCC, GGTCCC, 
AGTCCT, GGCCCT and AGCCTT belonged to the 
superior group for all 3 traits. 

These results indicate that g.6850+77 A>G (AA, AG), 
g.10213 T>C (TT, TC) and g.14047 C>T (CC, CT) 
provided the best SNP combination for CI 8:1, MUFAs and 



MS in Korea native cattle. The MDR method is a fast and 
reliable technique for multiple SNP analyses of potential 
applications in marker-assisted selection. In addition, this 
technique, a non-parametric statistical method for detecting 
gene interactions, is useful for small samples and performs 
an exhaustive search of all n-locus models, collapsing 
multi-locus genotypes into high- and low-risk classes 
(Ritchie et al., 2001). This suggests that, after further 
development, the MDR method may be used for the genetic 
assessment of quantitative traits. Finally, these genotypes of 
individual SNPs and their combinations may be important 
genetic markers for improving beef quality, which is one of 
the most important goals in the Korean cattle industry. 
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